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1. Introduction 

An interesting class of isotropic inhomogeneous optical media that are currently 
investigated in optics is represented by the so-called spherical media, i.e. the media with 
spherically symmetric distribution of refractive index. These media found application in 
theoretical proposals of many optical devices that provide remarkable optical properties 
within the framework of geometrical optics. For example, we mention the spherical 
gradient- index lenses embedded into an optically homogeneous region that provide 
perfect imaging of two concentric spherical surfaces. Their well known representatives 
are a Luneburg and an Eaton lens [H [2]- In addition, we mention the absolute 
instruments that provide perfect imaging of three-dimensional domains |3j. Their 
famous representative is a Maxwell's fish-eye |1] which, in fact, provides perfect imaging 
also in terms of wave optics 0. 

A common approach to designing the spherical medium of specific properties is 
based on solving an inverse scattering problem. This is considerably simplified due 
to the symmetry of the spherical media which enables to treat these media as two- 
dimensional, because each light ray lies in a plane through the center of symmetry. 
The classical methods for solving the inverse scattering problem have been derived for 
both the spherical gradient- index lenses as well as the absolute instruments [6]. 
Recently, an alternative method based on ideas of transformation optics has also been 
described [7]. This alternative method utilizes an equivalence of the central section 
of the spherical medium and a curved rotationally symmetric surface with constant 
refractive index referred to as a geodesic lens [H [9l [TOl [TT] . A great advantage of this 
method lies in a fact that the shape of a geodesic lens can be deduced relatively easily 
from the desired light ray trajectories for both the focusing spherical gradient-index 
lenses and the spherical absolute instruments. Then, once the shape of a geodesic lens 
is known, the refractive index of an equivalent spherical medium can be calculated by a 
straightforward procedure. 

In this paper we use this alternative approach based on the concept of geodesic 
lenses to design new optical devices referred to as multi-focal gradient-index lenses and 
multi-focal absolute instruments. These devices formed by isotropic spherical media 
with complex spatial structure provide different imaging properties for different light 
rays and, therefore, enable to achieve remarkable optical effects. In fact, we have already 
presented some multi-focal absolute instruments before [B] . However, now we show that 
the multi-focal absolute instruments as well as the multi-focal gradient-index lenses 
can easily be designed using the concept of geodesic lenses and that there is a close 
connection between them. As in our previous paper we perform the description 
within the framework of geometrical optics and we utilize the symmetry of spherical 
media to reduce the description to two dimensions. 

This paper is organized as follows. Section |2] is devoted to the derivation of the 
multi-focal gradient-index lenses embedded into an optically homogeneous region that 
provide perfect imaging of multiple pairs of concentric spherical surfaces. Section [3] is 
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aimed to the derivation of the multi-focal absolute instruments that provide different 
imaging properties for different spatial domains. Finally, section |4] concludes the article. 

2. Gradient-index lenses 

This section is devoted to the spherical gradient-index lenses embedded into an optically 
homogeneous region. We start with a brief summary of physics related to these optical 
devices and a procedure for solving the corresponding inverse scattering problem. Then 
we utilize the general results to design the multi-focal gradient-index lenses that provide 
perfect imaging of multiple pairs of concentric spherical surfaces. 

2.1. Gradient-index lens and equivalent geodesic lens 

Let us start with a brief summary of physics related to the spherical gradient-index 
lenses and equivalent geodesic lenses. 

Consider a plane with the system of polar coordinates (r, (f) that represents an 
optically homogeneous region with refractive index set to unity. Into this plane we 
place a gradient- index lens of unit radius with the center at the origin (see figure Q that 
is specified by a spherically symmetric refractive index n(r), which meets the boundary 
condition n(l) = 1. This gradient-index lens is, in terms of the transformation optics, 
equivalent to a geodesic lens which is a curved two-dimensional surface (waveguide) with 
rotational symmetry and refractive index set to unity (see figure [2]). We describe its 
shape by a radial coordinate p, an angular coordinate 6 and a value of function s{p), 
which represents the length of the surface measured along the meridian from the axis 
of symmetry to the given point. The coordinate transformation between the gradient- 
index lens and the geodesic lens follows from the comparison of equivalent optical path 
elements |7]. Assuming 9 = ip, we get 

p = nr, ds = n dr. (1) 

Clearly, the shape of a geodesic lens depends primarily on the form of function p(r). 
Within the section [2] we will assume that for r G [0,1] the function p(r) is increasing 
with a single maximum at r = 1 and, therefore, the shape of a geodesic lens is of the 
form sketched in figure |2} Although the function s(p) is unambiguous, it will be referred 
to as Si (p) to maintain consistency with section [s] 

Each light ray propagating in the spherical medium or on the geodesic lens 
is specified by a constant quantity L that is analogous to an angular momentum 
known from classical mechanics (therefore, we will call the quantity L simply angular 
momentum). It is given by 

L = nr sin a = p sin a (2) 

where a is an angle between the tangent to the ray trajectory and the radius vector in 
the spherical medium or, equivalently, an angle between the tangent to the ray trajectory 
and the meridian on a geodesic lens. The point of the ray trajectory where a = 7t/2 
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Figure 1: Geometrical configuration of tlie gradient-index lens embedded into an 
optically homogeneous region. 
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Figure 2: Geometrical configuration of the geodesic lens equivalent to the gradient-index 
lens. 



will be referred to as a turning point, its radial coordinate in the spherical medium as 
r+(L). Considering an element of the ray trajectory and denoting s[{p) = dsi(p)/dp, 
we obtain with help of ([T]) and ^ a differential equation 

Ldr Ls[{p)dp 



dip 



(3) 



that governs propagation of the light ray with angular momentum L in the spherical 
medium and on the geodesic lens, respectively. These formulas represent the necessary 
basis for mathematical formulation of an inverse scattering problem discussed in the 
following section. 



2.2. Inverse scattering problem 

The classical formulation of an inverse scattering problem for a gradient-index lens 
embedded into an optically homogeneous region lies in derivation of the refractive index 
n{r) from a given dependence of a scattering angle x{L) on the angular momentum L, 
where the scattering angle corresponds to the angle by which the incoming light ray is 
deflected from its original direction (see figure [T]). Here, we specify bending of the light 
rays in an equivalent, but slightly different way that is more suitable for addressing the 
inverse scattering problem related to designing the optical devices that provide stigmatic 
imaging of given surfaces. For a given light ray with angular momentum L, we choose 
a 'source' point Ps[rs(L), (y9s(-Z^)] on the straight line along which the light ray comes 
from infinity and we choose an 'image' point Pi[ri(L), y9i(L)] on the straight line along 
which the light ray continues to infinity after passing through the lens. In addition, we 
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introduce an angle A{p{L) which corresponds to an absolute value of the polar angle 
swept by the light ray during its propagation between the points Pg and Pi (see figure 
[T|. The notation suggests that the radial coordinates rs(-L), ri(L) as well as the angle 
Aip{L) are, in general, functions of L; this is necessary for considering the multi-focal 
lenses discussed later. These three functions fully determine bending of the light rays 
passing through the lens and, therefore, the inverse problem lies in derivation of the 
refractive index n{r) from the given functions rs(L), ri(L) and Aip{L). 

The mathematical formulation of this problem is based on expressing the polar angle 
A(f{L) and it results in an integral equation for the refractive index n(r) of the gradient- 
index lens or a similar integral equation for the function s[{p) that parameterizes a 
geodesic lens. The particular form of these equations depends on the form of function 
p{r) = nr. Under the assumptions for the function p{r) given above, these integral 
equations can be written in a compact form 



Ir+iL) ry/n'^r'^ - Jl p^ p^ - 

where the function g{L) is given by 

giL) = - ( A{p(L) + arcsin — — — + arcsin — 2 arcsin L ) . (5) 
2 V r,{L) n{L) ) 

The standard procedure for solving the integral equation for refractive index n(r) can 
be found in [1]. Here, we use an alternative approach that lies in solving the integral 
equation for the function s'^(p) and subsequent calculation of the refractive index n{r). 

The integral equation for the function s'i{p) is a kind of Abel integral equation and 
can be solved by a slight modification of the method presented in [B] . The general result 
that is independent of the particular form of function g{L) is 

, , , 2p d / fi g(L)>lL \ 

Once we known the function s^(p), we are able to calculate the refractive index n{r). 
Using the relations ([!]) we get 

r p 

Integrating this equation from the boundary of the lens to the turning point, given in 
the spherical medium by the radial coordinate r equal to r^{L) and on the geodesic lens 
by p equal to the angular momentum L, we obtain 



r+(L) = exp / s\{p)^\. (8) 




Since we have assumed that the function p{r) is monotonic, also the function r+(L) is 
monotonic and invertible. Therefore, in principle we can find the inverse function i^(r+) 
and use it for calculation of the refractive index by means of relation 

n(r+) = (9) 
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which holds for the turning point because there p = L. Omitting the lower index of the 
radial coordinate, we obtain the final result for the refractive index n{r) of the spherical 
medium forming the gradient-index lens. 

2.3. Multi-focal gradient-index lenses 

In the following we will describe the gradient-index lenses that provide stigmatic imaging 
of pairs of concentric spherical surfaces. 

In our previous paper |7j, we have discussed the simplest case when there is only 
one pair of these surfaces. This corresponds to the situation when all the light rays 
with L G [0, 1] that emerge from a point source located at a given surface and that pass 
through the lens are focused into an image point located at the second given surface. 
Such lenses can be found by solving the Luneburg inverse problem [1] and their well- 
known representatives are the Maxwell's fish-eye lens, the Luneburg lens and the Eaton 
lens. 

Here, we deal with the multi-focal gradient-index lenses that provide perfect 
imaging of multiple pairs of concentric spherical surfaces. For simplicity, we start with 
the case when there are two pairs of these surfaces, i.e. the light rays passing through 
the lens are divided into two bundles specified by subintervals of angular momentum 
[0, Li] and [Li, 1] which, in general, have the source and the image points located at the 
spherical surfaces of different radii. Then the radial coordinates rs{L) and ri{L) of the 
source and the image points, respectively, can be written as 

[ rs2 [ for L E (Li, IJ, 

where Vgi, rs2, Tii and are given constants greater than or equal to unity. Similarly, we 
write the angle Aip{L) swept by the light rays during the propagation from the source 
to the image in the form 

/ tor Le [O.Li] 

where Mi > M2 are positive constants (a more detailed analysis shows that for Mi < M2 
the integral equation Q has no solution). The different imaging properties of the 
light rays corresponding to two subintervals of L imply that the refractive index of the 
lens must have different functional dependence in two different regions separated by a 
boundary of radius r^{Li). Therefore, we denote 

^ ^ I ni2(r) for rG K(Li),l]. ^ ^ 

Similarly, a geodesic lens equivalent to the gradient-index lens must be parameterized 
gradually from the top downwards by two different parts of the function Si{p) formally 
written as 

''^'^ - [ si2{p) for pG[L„l]. ^''^ 
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Equipped with the introduced notation, we proceed to solving the inverse scattering 
problem. Using the general formula ^ we get the results for the functions s[^{p) and 
(p) in a compact form 



Sl2{p) = M + 



B, 



(14) 
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(15) 



Clearly, the results for the function Si2(p) and the parameters A2 and B2 are formally 
identical to the results presented in [7j as a solution of the Luneburg inverse problem. 
The refractive index n{r) of the gradient-index lens can be calculated by a general 



procedure described in section |2.2[ However, an analytical result can be obtained only 
for the special cases when the radial coordinates of the source and the image points are 
delimited to only two possibilities - unity and infinity, and also assuming that rsi = rs2 
and Til = ri2. Then the parameter Ai = 0, the parameter A2 becomes independent of p 
and we obtain an implicit formula 

2/B2 _ n^l/B2(„_^\A2/B2-l , {„_^\2A2/B2 



2r^/^'^{ni2ry'^'"'- ^ + [ni2r 
for the function ni2{r) and a similar formula 

r = p^2-Bi-i?2 U2 _ 







1 



P' 



B2 



(16) 



(17) 



which specifies the function ?2ii(r) = p(r)/r. However, in most cases the last equation 
can be solved only numerically. 

The described solution of the inverse scattering problem for the gradient-index lens 
that provide perfect imaging of two pairs of concentric spherical surfaces can easily be 
generalized for the cases with multiple pairs of these surfaces. Then the light rays passing 
through the lens are divided into bundles specified by subintervals of angular 
momentum that are separated by — 1 real constants < Li < L2 < . . . < L^^i < 1. 
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In the j-th subinterval of angular momentum, the functions rs{L), r\{L) and ^{L) are 
given by the quantities r^j, r^j and M^-tt, respectively, which represent generalization of 
the above notation. The refractive index n(r) is specified by the function nij{r) within 
an annular region r G [r+(Lj_i), r+(Lj)] and the corresponding part of an equivalent 
geodesic lens is parameterized by the function Sij{p) or by its derivative s'^^^p). Using 
the general formula (|6|), we find the function s'ij{p) in a compact form 



while we define = 1 and the parameters and are given as follows. For 
k G {1, 2, ... — 1}, the parameters Ak and Bk can be obtained from Ai and Bi listed 
in (15) by substituting of all indices 1 with k and 2 with k + 1. Similarly, for k = N 
we get the parameters Aj^ and B^ from A2 and B2 listed in (15) by substituting of all 
indices 2 with N. Finally, once the shape of a geodesic lens is known, the refractive 



index n(r) can be calculated by the procedure described in section 2.2 

The general formula (18) for the function s'ij{p) together with the formulas for the 
parameters A^ and Bj. obtained by generalization of (15) represent the main result of 
this section. It proves that the approach based on the concept of geodesic lenses is 
very efficient for designing the multi-focal gradient-index lenses - in any general case 
the shape of a geodesic lens can easily be deduced from the desired imaging properties 
and the refractive index of the gradient-index lens can consequently be calculated by a 
straightforward procedure. 



2.4- Examples 

Let us now present in figure [3] the ray tracing and the equivalent geodesic lenses for 
several examples of the multi-focal gradient-index lenses that provide perfect imaging 
of two concentric spherical surfaces. We focus on the special cases when the source and 
the image points are located at infinity or on the unit circle and we refer to the lenses 
according to the effective behavior of the light rays. 

As the first example we present a lens determined by parameters A\ = A2 = 0, 
i?i = i?2 = 1 which consists of a classical Maxwell's fish-eye in the outer region, but 
acts as a generalized Maxwell's fish-eye for the light rays passing through the inner 
region. The second example is a lens given by numbers = 0, = -B2 = ^, -Bi = 1 
corresponding to the combination of a classical and a so-called generalized Luneburg 
lens which focuses the light rays coming from infinity to two opposite points on the 
boundary of the lens. As the third example we present a lens with parameters A\ = 0, 
A2 = Bi = B2 = 1 which acts as the combination of an Eaton lens and an invisible 
sphere. 

The last example shown in figure |3] is probably the most interesting one. This 
lens effectively acts as the combination of a Luneburg lens and an invisible sphere and, 
therefore, provides an interesting optical effect. When looking through the lens, the 
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central region shows an image of the scene behind while the peripheral region shows a 
deformed image of the surrounding scene. This effect is apparent from the visualization 
shown in figure |4j 




(a) (b) (c) (d) 

Figure 3: Ray tracing and equivalent geodesic lenses for several multi-focal gradient- 
index lenses that effectively act as the combination of (a) a classical and a generalized 
Maxwell's fish-eye, (b) a classical and a generalized Luneburg lens (c) an Eaton lens and 
an invisible sphere (d) a Luneburg lens and an invisible sphere. The spherical medium 
is shown in light blue, the dashed circle has the radius r+(Li). 




Figure 4: A visualization of the optical effect provided by the lens shown in figure |3| 
that effectively acts as the combination of a Luneburg lens and an invisible sphere. 
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3. Absolute instruments 

Another interesting class of optical devices formed by the spherical media are commonly 
known as absolute instruments. They are characterized by providing perfect imaging 
of three-dimensional domains within the framework of geometrical optics which can be 
achieved when all the light rays propagate along closed trajectories. In this section, we 
briefly review the physics of the spherical media with spatially conflned light rays and 
we solve a corresponding inverse scattering problem. Then we use the general results 
to design the so-called multi-focal absolute instruments that provide different imaging 
properties for the light rays with different angular momentum. 

3.1. Spherical medium with spatially confined light rays 

Consider a plane with the system of polar coordinates (r, ip) and a spherical medium 
of radius R > 1 with the center at the origin that is specified by refractive index n(r). 
To ensure that all the light rays are spatially confined within this medium (see figure 
[5]) we assume that the corresponding function p(r) = nr has a single global maximum 
p(l) = 1, the value p(0) = and lirrir^R p{r) = 0. Clearly, the inverse function r(p) is 
multivalued; we denote the corresponding branches by r±{p). 

The spherical medium defined above is equivalent to a geodesic lens of the shape 
sketched in figure [6] Regarding the assumed form of function p(r), the geodesic lens 
must be parameterized by two branches of the function s{p) as follows. The branch Si{p) 
parameterizes the shape of its upper part gradually from the point Ti to the circle p = 1 
referred to as an equator. The branch S2(p) parameterizes its lower part gradually from 
the equator to the point T2, while we assume S2(l) = To maintain consistency 

with section |2] we assume that the upper and lower part of geodesic lens correspond in 
the spherical medium to the region inside and outside the unit circle, respectively. 

A light ray with angular momentum L propagates in this spherical medium within 
an annular region bounded by the circles r^{L) and r^{L) on which the turning points 
are located. Equivalently, it propagates on the geodesic lens within a region bounded by 
circles of radius p = L one above and one below the equator on which the turning points 
are situated. We denote an angle between two consecutive turning points by Aiptp{L). 
As being justified in [7], this angle remains the same for infinitely many geodesic lenses 
that can be mutually reshaped into each other and that are parameterized by functions 

s'lip) = <(p) + 

4(P) = <(p) - 4(P) (19) 

where the functions s'^{p) and s'^{p) describe the symmetric and the antisymmetric part 
of a geodesic lens, respectively. This represents a key fact that allows solving an inverse 
scattering problem discussed in the following section. 
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Figure 5: Geometrical configuration of tlie splierical medium witli spatially confined 
light rays. 




Figure 6: Geometrical configuration of the geodesic lens equivalent to the spherical 
medium with spatially confined light rays. 



3.2. Inverse scattering problem 

An inverse scattering problem for spatially confined light rays lies in derivation of the 
refractive index n{r) from the given function A(ftp{L). Expressing this angle by means 
of equation (|3]), we get integral equations 

r+{L) Ldr _ L {s[{p) - s'^{p))dp 

for the refractive index n(r) and the functions s\{^p^ and S2(p) that parameterize the 
geodesic lens. The procedure for solving the integral equation for the refractive index 
n(r) can be found in 0. Here, we use an alternative approach already used in section |2] 
that lies in solving the integral equation for the functions s'i(p) and S2(p) subsequent 
calculation of the refractive index nir). 



L{s[{p)-s'^{p))dp 



The procedure for solving the integral equation (20) starts by employing the 



relations (19). Then the integral equation is reduced to 
^1 L4(p)dp 1 



L pVp^ - 2 



Aip,,{L) (21) 
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where the only unknown function is s'^{p). This is again a kind of Abel integral equation 
and its general solution is 

This formula allows to calculate the shape of a symmetrical geodesic lens from a given 
function A(ptp{L). Then, by choosing the function s'^{p) that makes the geodesic lens 
asymmetrical, it is possible to design a variety of geodesic lenses that represent spherical 
media with a common function Aiptp{L). 

The refractive index n(r) of the spherical medium can be calculated by means of the 
method described in section [2l2j In fact, the described method is directly applicable to 



derive the refractive index within the unit circle. The refractive index outside the unit 
circle can be calculated by formally identical relations obtained by a simple substitution 
s'lip) ^ 4(P) and r+(L) ^ r_(L). 

3.3. Multi-focal absolute instruments 

Let us now utilize the general results for discussion of the absolute instruments that 
provide perfect imaging of three-dimensional domains. 

In our previous paper [7], we have discussed the simplest case when the angle 
Aiptp{L) swept by the hght ray during the propagation between two consecutive turning 
points is independent of L and equals to a rational fraction of vr. Then the ray 
trajectories are closed and each source point is sharply imaged into a set of image 
points depending on the character of ray trajectories. 

Here, we deal with the so-called multi-focal absolute instruments that provide 
different imaging properties for the light rays with different angular momentum. 
Analogously to section [2| we start with the case when the angle Ay9tp(-^) is piecewise 
constant on two sub intervals of L. We denote 

^^^^^'^^ = 1 ^.vr for (^3) 
where Bi > B2 are positive rational numbers and < Li < 1 (the condition for 



Bi and B2 is analogous to the condition Mi > M2 mentioned in section 2.3). This 
assumption implies that the refractive index of the spherical medium has different 
functional dependence within four annular regions separated by the boundaries of radii 
r+(L), 1 and r-{L), respectively. For the refractive index inside the unit circle we use 



the notation (12), in addition we denote 



n(r) = l for r G [l,r_(Li)] 

^ ' \ n2i{r) for r G [r^{L,),R]. ^ ' 

An equivalent geodesic lens is parameterized by two different parts of functions si{p) 



and S2(p). The former is already specified by (13), the latter is given by 



, ^^^alp) for P e [Li, 1] 

''^'^ ~ \ S2,ip) for pG[0,Li], ^2') 
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where the functions S22(p) and S2i(p) parameterize the lower part of a geodesic lens 
gradually from the equator downwards. Similarly, we split the functions Ss{p) and Sa(p) 
as 

^.w = I . i°' ""["/^jj (26) 

[ Ss2(p) [ Sa2(p) for p G [Li,l\. 

Equipped with the notation, we proceed to solving the inverse scattering problem. 



Using the general formula (22) we get the functions 



n2 



's2(p) = -^rH (27) 



'1 - {p/LiY VI -p 
B2 

that describe a geodesic lens symmetrical with respect to the equatorial plane. Choosing 
the functions s'^^i^p) and s'.^2i.P) '^^ S^t a variety of asymmetrical geodesic lenses that are 
equivalent to the multi-focal absolute instruments. 

In analogy with the results of section |2.3[ we focus on a class of geodesic lenses 



given by the choice s'^^(p) = Ai and s^2(p) = -^2, where Ai and Ai are non-negative 
real numbers. Then the geodesic lens is described by a set of functions 

B2 



s'uip) 


= ^1 + 


Bi 


-B2 










ip/Liy 




= A2 + 


B2 




Vl- 


p' 


42 (p) 


= A2- 


B2 






p' 


4i(p) 


= A- 


Bi 


-B2 











+ 



2 



Vl -p^ 

In such case, the refractive index of the spherical medium can be calculated by means of 



equations derived in section 2.3 The functions ni2{r) and n22{r) are given by a formula 



obtained from (16) by a simple substitution A2 — ?■ A2 and B2 — ±B2- Similarly, the 



functions nii(r) and n2i(r) are given by an equation obtained from (17) by a substitution 
Ai Ai, Bi ±{Bi - B2) and B2 ±62- 

Similarly as in section [2| the described solution of the inverse scattering problem 
for the absolute instruments that provide different imaging for the light rays specified 
by two subintervals of angular momentum can easily be generalized for the cases with 
of these subintervals separated by — 1 real constants < Li < L2 < . . . < L^-i < 1. 
Then the function A{ptp{L) is given in the j-th subinterval of angular momentum by 
the quantity Bjti, which represents generalization of the above notation. The refractive 
index is given by the functions nij(r) and n2j{r) inside and outside the unit circle, 
respectively. An equivalent geodesic lens is parameterized by the functions sij(p) and 
S2j{p) above and below the equator, respectively, and its symmetric and antisymmetric 
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Ss,ip) = T. / . (29) 



part are described by functions Ssj{p) and Saj(p), respectively. Using the general formula 
(22), we find the function s'gj{p) in a compact form 

k=j ^1 - {p/Lk 

while we define L^v = 1 and -Bat+i = 0. Choosing the functions s'^j{p) we get a variety 
of asymmetrical geodesic lenses and, finally, we calculate the refractive index of the 
corresponding multi-focal absolute instruments by the procedure described in section 

The general formula (29) represents the main result of this section. It proves that 
the approach based on the concept of geodesic lenses is very efficient also for designing 
the multi-focal absolute instruments. Moreover, similarity of the general formulas ( 18 ) 
and (29) reveals that there is a close connection between the multi- focal gradient-index 
lenses and the multi-focal absolute instruments. 



3.4- Examples 

Let us now show in figure [7] the ray tracing and the geodesic lenses corresponding to 
several examples of the multi-focal absolute instruments that provide different imaging 
properties for the light rays specified by two subintervals of angular momentum. 

As the first example shown in figure [7||^a) we present an absolute instrument specified 
by parameters Ai = A2 = 0, Bi = 2 and B2 = I- This corresponds to the simple 
case when the geodesic lens is symmetrical with respect to the equatorial plane. Its 
central part corresponding to p > Li is formed by a spherical surface on which the light 
rays propagate along great circles, the parts corresponding to p < Li one above and 
one below the equator are formed by more complex surfaces on which the light rays 
propagate along the geodesic curves that cannot be simply classified. The refractive 
index of an equivalent spherical medium is given by the profile of Maxwell's fish-eye 
within an annular region bounded by the dashed circles of radii r±{Li) and by the 
numerically calculated profiles outside this annular region. The light rays emerging from 
a point P are divided into two color-coded bundles corresponding to two subintervals of 
L. The red light rays with L G {Li, 1] propagate along the circular trajectories within 
the dashed bounded annulus and they meet at the image point Q before they reach the 
starting point P again. On the other hand, the blue light rays with L G [0, Li] propagate 
also outside the dashed bounded annulus and they meet at the image point T before 
returning to the point P. The blue and the red trajectories also differ in the polar angle 
between two consecutive turning points according to the choice of parameters Bi and 
B2, respectively. 

The second example shown in figure [Tj^b) is specified by parameters ^1 = 0, ^2 = ^, 
Bi = 1 and = |. The geodesic lens is asymmetrical with respect to the equatorial 
plane and the refractive index is given by a Luneburg profile within the dashed bounded 
annulus. The red light rays with L G {Li, 1] emerging from a point P propagate along 
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(a) (b) (c) 

Figure 7: Ray tracing and equivalent geodesic lenses for the multi-focal absolute 
instruments given by the following parameters: (a) Ai = A2 = 0, Bi = 2, B2 = 1, 
(b) Ai = 0, ^2 = Bi = 1, B2 = |, (c) Ai = A2 = B2 = 1, Bi = 2. The black 
circle corresponds to the circular trajectory at r = 1 with maximum possible angular 
momentum L = 1, the dashed circles have the radii r±[Li). 

the concentric ellipses and meet at the image point Q before returning to the point P. 
The blue light rays with L G [0, Li] emerging from the point P meet at the image point 
T that is different from Q and the angle between the consecutive turning points is twice 
that of the red rays. 

Finally, we show in figure the third example given by parameters Ai = A2 = 
B2 = 1 and Bi = 2. In this case, the refractive index is given by an Eaton profile within 
the dashed bounded annulus and, therefore, the light rays propagate within this region 
along the confocal ellipses. The red rays with L G (-^^1,1] emerging from a point P 
return to this point after making one loop around the origin. On the other hand, the 
blue rays with L G [0,Li] emerging from the point P need to make two loops around 
the origin to close their trajectory. 
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4. Conclusion 

In this paper we have presented a general approach to designing isotropic spherical media 
with complex spatial structure that provide different types of imaging for different light 
rays. This approach based on equivalence of the spherical medium and the corresponding 
geodesic lens proved to be very efficient for designing the multi-focal gradient-index 
lenses embedded into an optically homogeneous region as well as the multi-focal absolute 
instruments that provide perfect imaging of three-dimensional domains. 
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